Chaotic advection and targeted mixing 
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The advection of passive tracers in an oscillating vortex chain is investigated. It is shown that 
by adding a suitable perturbation to the ideal flow, the induced chaotic advection exhibits two 
remarkable properties compared with a generic perturbation : Particles remain trapped within a 
specific domain bounded by two oscillating barriers (suppression of chaotic transport along the 
channel), and the stochastic sea seems to cover the whole domain (enhancement of mixing within 
the rolls). 



PACS numbers: 05.45.Ac, 05.60. Cd, 47.27.Te 

In order to enhance mixing in flows one often relies on 
chaotic advection Q. This phenomenon translates the 
fact that despite the laminar character of a given flow, 
trajectories of fluid particles or advected passive tracers 
are chaotic. As a consequence, mixing is considerably 
enhanced in regions where chaos is at play, as it does 
not need to rely on molecular diffusion. Such phenom- 
ena are observed in a wide range of physical systems j3| 
and have fundamental applications, for instance, in geo- 
physical flows 0, Q| or chemical engineering. When con- 
sidering two dimensional incompressible flows, chaotic 
advection occurs when the flow is unsteady, see for in- 
stance Refs. 0, IE 0, 0|- One pecuharity of these flows 
is that the dynamics of passive particles can be tack- 
led from the Hamiltonian dynamics point of view. The 
canonical variables are directly the space variables, hence 
the phase space is formally the two dimensional physical 
space. This feature allows a direct visuahzation of the 
phase space in experiments, and in a certain sense makes 
this framework an ideal one in order to confront theoret- 
ical results on Hamiltonian dynamics with experiments. 

In this Letter, we address the problem of targeted 
mixing. Namely it is often desirable to reduce chaotic 
transport (see Ref. 0) while one also often wants more 
chaos in order to enhance mixing (see, for instance, 
Refs. for microfluidic and microchannel devices). 

In order to achieve such properties, we consider the dy- 
namics of passive tracers in an array of alternating vor- 
tices. This flow has been realized experimentally using 
Rayleigh-Benard convection or in a more controlled way 
by using electromagnetic forces 0, 0, 0|- The primary 
interest in this flow resides in the fact that being gener- 
ated by quite a few hydrodynamic instabilities, it may 
be considered as one of the founding bricks of turbu- 
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lence. As such, understanding its influence on the advec- 
tion of passive or active quantities is considered a nec- 
essary flrst step in order to uncover the different mecha- 
nisms governing for instance front propagation in turbu- 
lent flows EUll. 

We consider the following integrable stream function 
which models the experiment with slip boundary condi- 
tions : 



'^i){x,y) = a sin a; sin y, 



(1) 



where the x-direction is the horizontal one along the 
channel and the y-direction is the bounded vertical one. 
The constant a is the maximal value of vertical velocity. 
Passive particles follow the streamlines of 5*0 depicted on 
Fig. ^a). The dynamics given by Eq. (QJ is integrable. 
Therefore no mixing occurs. The fluid is limited by two 
invariant surfaces y = tt and y = corresponding to the 
top and bottom roll boundaries. The flow has hyperbolic 
flxed points on these two invariant surfaces which are lo- 
calized at a; = rmr for to G Z. These points are joined by 
vertical heterochnic connections for which the stable and 
unstable manifolds coincide. The phase space, which is 
here the real space, is then characterized by a chain of 
rolls with separatrices localized at a; = rmr for m e Z. 

In the experiments [5, ^ a typical perturbation 
f{x,y,t) is introduced as a time dependent forcing in or- 
der to trigger chaotic advection and then to study the re- 
sulting transport and mixing properties. More precisely 
the perturbation modifles the stream function as : 



■fc{x,y,t) = ^o{x + f{x,y,t),y). 



(2) 



For instance, the following stream function has been pro- 
posed to model an experimental situation : 



'^i{x, y,t) — a sm{x + e sin wpi) siny. 



(3) 



i.e. the perturbation is / = esiwujQt and describes the 
lateral oscillations of the roll patterns where e and ujq are 
respectively the amplitude and the angular frequency of 
the lateral oscillations. Without loss of generality, we 
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Figure 1: Poincare sections (a) for the stream function JJl 
which are also the streamlines and (&) for the stream func- 
tion lOl. The parameters are a — 0.6 and e — 0.63. 



assume that ujq — I. The resulting streamhnes which 
are again closed curves correspond to lateral oscillations 
in the x-direction of the streamHnes depicted in Fig.nia) 
with a periodic displacement of — esini. However chaotic 
advection is triggered. 

In order to infer transport and mixing properties we 
use Poincare sections of some trajectories of tracers. A 
typical Poincare section associated with the stream func- 
tion 5'i given by Eq. JSJ is represented in Fig. Wib) for 
e — 0.63 and a = 0.6. It shows that the transport along 
the channel is greatly enhanced @ compared with no 
transport in Fig. Qla). The advected tracers have now 
a non integrable Hamiltonian dynamics. The reason be- 
ing that, under the perturbation, the vertical heteroclinic 
connections between the vortices break down and the 
stable and unstable manifolds intersect transversely thus 
generating chaotic advection of passive tracers along the 
channel [lj|. Therefore the mixing properties of the flow 
are enhanced. However, most of the time islands of sta- 
bility remain, among which the ones located around the 
centers of the vortices. They are characterized by regular 
(quasiperiodic) trajectories and a chaotic region around 
and between the rolls. Mixing within the center of the 
rolls can only rely on molecular diffusion while tracers in 
the chaotic sea are advected unboundedly along the chan- 
nel and chaotic mixing occurs. The main observation is 
that the regular patterns persist and chaos or mixing are 
not well developed. 

In order to obtain what we called targeted mixing, we 
propose to adjust the time-dependent forcing / which we 
assume to depend on y and t such that the stream func- 



tion ^'c given by Eq. ^ has the following properties : 
First its attendant chaotic transport along the channel 
(a;-direction) is confined by impenetrable barriers. Sec- 
ond, when restricted to each of these independent cells, 
we want to obtain good mixing properties. 

A first possibility to create barriers given by invari- 
ant tori of the dynamics would be to use the method 
proposed in Refs. for controUing Hamiltonian sys- 

tems. This method uses an affine control (modifying the 
Hamiltonian or the stream function additively) . However 
applying this method to the system under consideration 
would create some barriers but would also regularize the 
inside of the cell, thus drastically reducing the mixing. 

In order to proceed in finding a more suitable pertur- 
bation f{y,t), we start from a more generic situation for 
which the stream function given by Eq. 1^ is a particu- 
lar case. We consider a generic integrable stream func- 
tion '^o{x,y) which is L-periodic in x, and such that 
d^'^o{x,0) = d^'^oix,H) = for all 2: € R where H 
is the height of the channel. The perturbed stream 
function is written as in Eq. ^ where f{y,t) is such 
that (i) there are barriers aX x = Xk{y,t) (for k £ 2.) 
which prevent the chaotic advection along the x-direction 
and (m) there is a high mixing inside a region bounded 
by two consecutive barriers, e.g., xi and X2. We no- 
tice that the perturbed stream function is obtained by 
adding the perturbation to the canonical variables of the 
stream function and not to the stream function itself as 
m Refs. 0,0. The equations of the barriers are written 
as a; = Xk{y,t) = kL + ip{y,t) for /c g Z where if is a 
function to be determined. Using Hamilton's equations 

X ^ ^ -^^{x + f,y) - ^^P^ix + f,y) , 



dy dy 



dy dx 



and by imposing that Xk locates a barrier we get for the 
particles on the barrier 

df . dip 
dt y dy 

Thus a possible solution is obtained when f + ip is only 
a function of time, denoted and the condition on ip 
is 



dip _ d^o 
dt dy 



The above equation is solved as 

ip>{y,t) = -rdy-^omt),y), 

where the linear operator F is a pseudo-inverse of dt, i.e. 
acting on v{y,t) = Y^k'^k^'^* 



ik 



v,^,kt_ 



fc#0 



3 




Figure 2: Streamlines of the stream function lO at (a) t = 
and (6) t — 37r/4. The parameters are a — 0.6 and e = 0.63. 

The perturbation / is then 

f{y,t)^m+^dy'fo{m,y), (4) 

where <& is any function of time. 

We now apply these results to the array of vortices 
given by the stream function In order to remain as 
close as possible to the experimental setup, we consider 
= esint. The modified stream function given by 
Eq. 1^1 is then 

"i'dx, y,t) = a s\n[x + e sinf + a cos yCc{t)] sin j/, (5) 

where 

^^(*) ^ ^;^^2n+i(e) cos(2n + 1)^, (6) 

n>0 " 

and Ji (for I e N) are Bessel functions of the first kind. 
We notice that y — Q and y — n are still boundaries for 
the streamfunction 10 which comes from the fact that 
the modification is only on the x-term of the original 
stream function. The streamlines are then slightly mod- 
ified (non-uniformly in y) as it can be seen on Fig. I^a) 
and (6) which depict the streamlines of the stream func- 
tion © for e = 0.63 and a = 0.6 at two different times 
t = and t = 37r/4 respectively. Moreover, the displace- 
ment of these rolls remains parallel to the x-axis as it is 
for the stream function given by Eq. (jSJ. The fact 
that the streamlines of the stream function ^ and the 
ones of Eq. @ look similar comes from the fact that the 
stream function '^c is a small modification of (for a 
small) since we have |5'c(a;, y, t) — ^i{x, y,t)\ < 




Figure 3: Poincare section for the stream function JSJ. The 
parameters are a = 0.6 and e = 0.63. 

Nevertheless, the dynamics of tracers is completely dif- 
ferent in both cases and targeted mixing is achieved. For 
the same values of e and a as in Fig. 016), a Poincare 
section of the dynamics of passive tracers in a fiow de- 
scribed by the stream function ^'c given by Eq. Q is 
represented in Fig. We notice that there are invari- 
ant surfaces which have been created around x = (mod 
27r) (bold curves) and that mixing within these barriers is 
considerably enhanced. The equations of these transport 
barriers along the x-direction are : 

X — Xk{y,t) =^ 2kn — acoayCe{t). (7) 

To be more specific, the barriers Q are degenerate in- 
variant tori since each of them is a heteroclinic connexion 
between two hyperbolic periodic orbits (with period 27r), 
e.g., one locate at y = and x{t) = —aCe{t) and the 
other one at y = tt and x{t) — aCe(i), moving in oppo- 
site directions along the channel. Regarding the enhance- 
ment of mixing inside the cell bounded by the two barri- 
ers observed in Fig. El we notice that most of the regular 
trajectories observed for the stream function Vfi are bro- 
ken by the perturbation (in comparison with Fig.^6)). 
It means that, the perturbation given by Eq. Q cre- 
ates two invariant surfaces around x — and x — 27r, 
and destabilizes the inside of the cell and in particular 
the remainders of the regular motion near x — tt/2 and 
X = 37r/2. 

In order to get some clues on the origin of mixing 
enhancement, we consider the region in the middle of 
the cell of length L in the x-direction and compute the 
Melnikov function which measures the distance between 
the stable and unstable manifolds. First, we perform 
the (time-dependent) canonical transformation {x,y) i-^ 
{x\y') generated by F2{y,x',t) = yx' -T'i/o{'^{t),y). In 
the case of given by Eq. 1^ , the stream function be- 
comes 

^c{x', y\t) ^ a[sin{x' + $(t)) - sin($(t))] sin y' . 
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Figure 4: Numerical simulation of the dynamics of a dye at 
t = 30, t ^ 50, t ^ 70 and t = 140 (from top to bottom) : 
left column for the stream function 101 and right column for 
the stream function ((SJ. The parameters are a — 0.6 and 
e = 0.63. 

The Melnikov function is defined as Q 

/ + (X> 
{*o,*J(xo(t-to),z/o(i "to))d<, 
-oo 

where {xo{t),yo{t)) is a solution of the unperturbed sys- 
tem ^fQ. The Poisson bracket between and is given 
by 

{^'o, ^c} = 0^(1 — COS a;) sin($(f)) sinj/cosy. 

It is straightforward to see that the Melnikov function 
vanishes when x — mod 2tt, i.e. on the boundaries of 
the cell, and is maximum (in absolute value) when x — 
TT mod 27r, i.e. in the middle of the cell. Therefore it is 
expected that the mixing is enhanced since then, there is 
a maximum fiux inside the cell 

In order to test the robustness of the method and to 
try an experimentally more tractable perturbation, we 
truncate the series given by Eq. The simplest possi- 
ble perturbation is obtained by considering only the first 
term of the series C^it) in Eq. © which leads to the 
following stream function 

y,t) = a sin (x + e sint — 2aJ^i{e) cos y cost) sin y, 

(8) 

which has a much simpler time-dependence. We per- 
formed a Poincare section for this simplified stream func- 
tion given by Eq. ijHJ for a — 0.6 and e — 0.63. 
This Poincare section looks identical to the one depicted 
on 13 showing that the simplified perturbation remains 
efficient. 



In Fig. 31 we depict a numerical simulation of the dy- 
namics of a dye in the fiuid. The left column shows the 
evolution of the tracers for the stream function ^'i given 
by Eq. lO . The right column shows the mixing of a dye 
within a cell which is limited by two barriers created by 
the stream function '^c given by Eq. We see that the 
mixing occurs through a combination of stretching and 
folding of the dye. The corresponding figure for '^s given 
by Eq. JHJ shows some advected particles outside the cell 
but the same picture inside. 

In conclusion, we considered a situation of time- 
dependent oscillating vortex chain. We have shown by 
adding a suitable perturbation (forcing) it is possible to 
suppress the chaotic advection along the channel by cre- 
ating dynamical barriers among the vortex chain, and at 
the same time, to enhance the mixing of passive tracers 
inside each of these isolated cells. The effect is obtained 
in spite of the oscillations of the rolls and is not too much 
affected by the use of a simplified version of the exact 
stream function. 
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